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DIFFERENTIAL MANEUVERING SIMULATOR DATA 


REDUCTION AND ANALYSIS SOFTWARE 

By Gary P. Beasley and Richard S. Sigman* 

Langley Research Center 

SUMMARY 

A multielement data reduction and analysis software package has been developed 
for use with the Langley differential maneuvering simulator (DMS). This package, which 
has several independent elements, was developed to support all phases of DMS aircraft 
simulation studies with a variety of both graphical and tabular information. The overall 
software package is considered unique because of the number, diversity, and sophistica- 
tion of the element programs available for use in a single study. 

The initial use of the software package was in support of a joint NASA/Navy fighter 
study. The purpose of this paper is to discuss the overall DMS data reduction and analy- 
sis package, as used in the fighter study, by reviewing the various elements of the soft- 
ware, showing typical results obtained from the study, and discussing how each element 
was used to support the study. 


INTRODUCTION 

Langley Research Center (LRC) has in operation the six-degree-of-freedom differ- 
ential maneuvering simulator (DMS) (refs. 1 and 2). This simulator was designed to be 
used to evaluate the differential performance of current and proposed aircraft. In order 
to effectively use information being obtained from the DMS, it was necessary to develop 
a data reduction and analysis capability. To aid in determining what should be included 
in this capability, a number of areas of interest were reviewed and a number of groups 
were consulted. Those consulted included National Aeronautics and Space Administration 
(NASA) design teams, military and NASA pilots, energy management and weapons systems 
experts, and various other military and civilian groups. The data reduction and analysis 
requirements resulting from this consultation and review were developed by NASA and 
used as a guideline in the development of the necessary software elements. 

Prior to the development of a satisfactory data reduction and analysis package, a 
joint NASA/Navy study was initiated. In order to have the complete DMS data reduction 
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and analysis software available for this study in the shortest possible time, personnel 
from the Naval Ship Research and Development Center (NSRDC) were consulted and uti- 
lized. The NSRDC group had had considerable programing experience along the lines of 
the DMS requirements. This group prepared new programs and modified existing NASA 
programs to meet the extensive DMS data reduction and analysis requirements. Some 
additional modifications and software development were performed by LRC and Naval 
Air Development Center (NADC) personnel. All the programs were designed to utilize 
as input either DMS input data used to define the aircraft simulated or DMS data resulting 
from simulation. 

The initial use of the software package was in the joint NASA/Navy simulation 
study started in July 1971. The main purpose of the study was to evaluate the differential 
maneuvering performance of several fighter-type aircraft using air combat maneuvering 
(ACM) as the prime task. The airplane definitions were provided by the Navy and the 
Defense Intelligence Agency (DIA) and were programed on one of the Control Data 6600 
computer systems by LRC personnel. The simulated aircraft were flown by current 
Navy fighter pilots experienced in ACM. 

Included in this paper are a discussion of the various elements of the DMS data 
reduction and analysis software package, a presentation of typical output provided by the 
different elements, and a brief review of how each element can be used. The information 
presented in this paper is not intended to be a user's guide but is an introduction to each 
program. 


SYMBOLS 


Cj) drag coefficient 

Cl lift coefficient 

E s specific energy 

g acceleration due to earth gravity 

h altitude (h = -Z) 

M Mach number 

n load factor 

probability of kill 
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specific excess power 


P,Q,r 

R 

t 

V 

X,Y,Z 

AX,AY,AZ 

x,y,z 

a 

J8 

e 

A 

^e 

X u/D 

A W 

a H 

d V 


roll, pitch, and yaw angular velocity, respectively 


range between aircraft 


time into run 


total aircraft velocity 


coordinates of inertial axis system 


difference in inertial coordinates, respectively, of two aircraft 
coordinates of body axis system 
angle of attack 
sideslip angle 

deviation angle, cos-1 XAX + YAY + ZAZ 

|R|fv| 

off^boresight angle (cos~l(x/R)) measured in body axis system, 

lAa - 4 = ft; I /d + x w 

component of X in x-y plane, tan - l(y/x) 

component of X with respect to x-y plane, sin _ l -p=^== 

\jx 2 + y2 

component of X in x-z plane, tan _ l(-z/x) 
component of X with respect to x-z plane, tan - * — 


(/z 2 + x 2 


angular rate of change of range vector in X-Y plane 


angular rate of change of range vector with respect to X-Y plane 



total angular rate of change of range vector in inertial coordinate system, 


a 


T 



<p,9,i(/ Euler roll, pitch, and yaw angle, respectively 


Subscripts: 

A,B aircraft located in DMS sphere A and B, respectively 

x,y,z along x-, y-, and z-axis, respectively 

max maximum 

Dot over a symbol denotes derivative with respect to time. 

STUDY BACKGROUND 

The NASA/Navy simulation study used the Langley differential maneuvering simu- 
lator (DMS) shown in figure 1. The DMS has a representative instrumented cockpit in 
each sphere and each cockpit has its own throttle, control stick, buffet system, and other 
instruments. The pilot in his cockpit is presented a projected visual scene, which repre- 
sents the sky, earth, and horizon, and a target which simulates the other aircraft. The 
projected target moves relative to the sky earth scene in response to each pilot's control 
inputs. In addition, the pilots are given instrument readings (attitude, altitude, airspeed, 
etc.) indicating his aircraft's status. Standard aircraft controls are used to provide inputs 
to one of the Langley CDC 6600 computer systems. The computer system performs all 
computations by using the aircraft descriptions programed and provides output which 
drives the DMS projectors and instrumentation. All of this is done in a real-time mode. 
More detailed discussions of the DMS are given in references 1 and 2. Figure 2 shows 
schematically the major elements involved in the simulation. 

Figure 2 also indicates the data output channels used by the DMS during a simula- 
tion. The major output of the DMS during the NASA/Navy study was the DMS Data Tape. 
Fifty-six variables are sampled and stored on a disk at a preselected rate (every 0.5 sec 
during the NASA/Navy study) during each computer run. At the end of a series of runs, 
the information on the disk was transferred to magnetic tape which was used as a perma- 
nent record for further data reduction and analysis. The variables sampled and stored 
on the DMS Data Tape are shown in table I. Other output channels used during the study 
included eight channel strip chart recorders, which monitored selected variables during 
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a run, and a cathode ray tube (CRT) display giving real-time three-dimensional trajectory 
plots of each aircraft. 


SOFTWARE REQUIREMENTS 

The data reduction and analysis requirements of the DMS can generally be assigned 
to one of three different categories. These categories include (1) pretest programs 
which provide information for pilot briefings (expected trends and areas of potential air- 
craft superiority or inferiority), (2) daily used data reduction programs which provide 
information on individual runs prior to the next day's runs, and (3) posttest programs 
which provide analysis of sets of runs for aircraft combinations as a function of pilots, 
initial condition (I.C.), and other parameters. The three types used in support of the 
DMS and the major elements of each are briefly reviewed in the following sections. The 
review consists of a discussion of the purpose and main components of the individual pro- 
grams and a description of the data input required and the output generated with examples 
of typical output. The information presented is not intended to be a user's guide to the 
programs but is instead an introduction to each program. Detailed information on pro- 
gram operation must be obtained from program operators. 

| 

Pretest Programs 

’ ' i 

Two particular types of information are provided by pretest programs for use inde- 
pendently or in support of DMS studies. The first type includes programs which provide 
digital simulation of ACM to evaluate possible maneuvers and to establish trends in offen- 
sive and defensive times for different aircraft combinations. The second type considers 
specific excess power and provides a theoretical indication of relative energy rate advan- 
tages and disadvantages. Both types of programs are discussed in this paper. 

Digital ACM simulation .- The two digital ACM programs which are available to 
support DMS studies were both designed to accomplish several main objectives. These 
programs were not used in the NASA/Navy study but can be used to conduct aircraft and 
weapons systems design trade-off studies and to investigate maneuvers for combinations 
of new aircraft-weapons systems. 

Analytical simulation of tactical air combat (ASTAC) .- Program ASTAC, which is 
a modified version of the TAC Avenger computer program (ref. 3), provides a compari- 
son of two aircraft in hard ACM situations and indicates the outcome of the maneuvering 
for a given I.C. and aircraft combination. Maneuvers used in the program are based on 
fighter pilot comments of what they would do under specific ACM situations on the basis 
of the two aircraft's relative positions, rates, etc. The major input data required for the 
ASTAC program includes the same drag, lift, thrust, and other aerodynamic data used in 
the DMS aircraft definitions. The ASTAC input also contains conditions and limits for 
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available maneuvers and various radar, gun, and missile characteristics for evaluating 
radar and weapons use and effectiveness during ACM. 

The output currently being obtained from ASTAC is summarized in table II. As 
indicated in the table the output includes single run data which (1) gives a time history of 
specific aircraft information, (2) indicates what maneuvers are selected, and when, for 
each aircraft, and (3) prints gun bursts and probabilities of kill for each gun burst. At 
the end of a run,, the program outputs the time on offense for each aircraft, information 
on the speed, load factors, and altitudes encountered during the run, and a summary of 
gunfire and missile probabilities of kill. In addition to the single run data, a multirun 
output is provided which summarizes the gun and missile information for all runs and 
provides average values of this information. 

Adaptive maneuvering logic (AML) .- Program AML was developed under contract 
to LRC by Decision Science, Incorporated (DSI). The primary objective of this program 
was to develop a method by which maneuvers for the various situations arising in air-to- 
air combat could be obtained without relying on pilot expertise or retrospective opinion; 
that is, the method would find suitable tactical maneuvers for any given situation of one's 
own and the opponent's aircraft by using only the basic capabilities of the two aircraft 
and weapons systems. 

Two characteristics distinguish Program AML from most other air-to-air combat 
simulation programs including ASTAC . The first feature of the program is its capability 
to adapt its maneuvers to changing characteristics of one's own and/or the opponent's air- 
craft. The second feature of the program is that it is not based on long-term classical 
maneuvers, such as "high-speed yo-yo" or "scissors" which are normally used to describe 
air-to-air combat. In contrast, the basic AML program maneuvers are elemental and 
utilize the following capabilities; 

(1) Change in bank angle 

(2) Change in load factor 

(3) Change in thrust 

The major inputs to the AML program include the thrust, drag, lift, maximum and 
minimum Mach numbers as a function of altitude, and maximum allowable load factor as 
a function of altitude and Mach number for both aircraft. This information is again the 
same as that used in the DMS aircraft definition. 

The operation of the AML program is as follows. The two aircraft start at a given 
initial condition and the maneuvers that follow are determined at decision points occurring 
at equal time intervals into the program. At each such decision point, the program maps 
the present physical situation of the two opposing aircraft into a given cell of a situation 
matrix. It then calculates a value associated with that particular cell. Next, it deter- 
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mines a number of elemental trial maneuvers which are likely candidates for successful 
maneuvers for that present situation. Then, by predicting the flight path of one aircraft 
for the trial maneuvers and by extrapolating the flight path of the opponent's aircraft for 
this time period, it determines the cells of the situation matrix in which the trial maneu- 
vers would end and the value of each of these cells. That trial maneuver which obtains 
the highest cell value is then executed. The questions that are answered in order to 
determine the cell value of the situation matrix are given in table EH. Each question is 
assigned a weighted score that is used to determine cell values. 

The types of output provided by the AML program are shown in table, IV. A more 
complete discussion of the AML program is presented in reference 4. 

Specific excess power programs .- Specific excess power P s of an aircraft is 
defined as the difference between the thrust available and the total drag multiplied by the 
aircraft velocity and divided by the aircraft weight. This gives a measure of the rate at 
which an aircraft is gaining or losing energy and the ability of an aircraft to change energy 
levels. This type of information, when used relative to another aircraft, indicates areas 
in which a maneuvering advantage exists. Specific excess power P s data for the 
NASA/Navy study were used to brief pilots on potential areas of superiority and inferi- 
ority for a particular aircraft combination prior to the actual simulated ACM. A sche- 
matic illustrating the P s type of programs available is included in figure 3. A brief 
description of each of these programs follows. 

Aircraft maneuvering performance (AMP): Program AMP was developed to gen- 
erate energy maneuverability data for a single aircraft by using the same lift, drag, 
thrust, and maximum and minimum altitude information as were used in the DMS. The 
data input to the program includes Cl as a function of Mach number M, angle of attack 
a, and altitude h; Cj) as a function of M and Cl; thrust as a function of M and 
h; Cl max as a function of M; and maximum and minimum altitude as a function of 
M. This information is in the same format as that used in the DMS. 

f 

The program operates by first selecting values of n z , M, and h as specified 
by the user. Angle of attack is iterated until the lift generated is equal to or greater 
than that necessary for the given n z , M, and h. Once this condition is met, the spe- 
cific excess power P g for that M, n z , and h is calculated. During the calculations, 
the program checks to determine whether Cl max or flight envelope limits have been 
exceeded. 

The output of the program consists of line printer output which lists P g and the 
iterated a at specified values of n z , M, and h. In addition to the line printer output, 
a magnetic tape is generated which contains the same data. The magnetic tape is saved 
for a permanent record and is used as input for a number of other P s programs as 
shown in figure 3. 
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Program CONTOUR: Program CONTOUR takes P s , h, and M data for a single 
aircraft and arranges the information for the plotting of a series of constant P s con- 
tours. The only input to the program is the tape generated by Program AMP which con- 
tains all the necessary information for a given aircraft. 

Program CONTOUR operates by taking the AMP tape for a particular aircraft and 
placing P s and a. information into arrays according to Mach number and altitude. 
These arrays are then used to find a set of P s contours for a specific load factor by 
selecting a Mach number and iterating altitude until a P s corresponding to one of the 
desired contours is found. At this point the h, M, and P s values are placed in arrays 
which are used to determine coordinates for plotting the contours or data for printout. 

This process is repeated for a series of Mach numbers. After a set of contours is 
obtained for a particular load factor, the data are stored on magnetic tape for later use, 
the load factor is changed, and a new set of contours is obtained. The number of load 
factors tested is a variable. The plot data stored on tape from CONTOUR are used as 
input t6 one of the Langley plotting systems. A typical plot resulting from CONTOUR is 
given in figure 4 . This figure shows a set of P s contours for a load factor of 1 plotted 
at various altitudes and Mach numbers. 

Program NEMO: Program NEMO creates a plot file of P s as a function of turn 
rate and turn radius at given altitudes and Mach numbers for two aircraft. This plot file 
can be processed to obtain plots on 30 or 16 inch (0.76 or 0.41 meter) paper. 

Input to the program consists of two tapes and two data cards. The tapes, contain- 
ing P g data, are created by Program AMP and each tape represents a single aircraft. 
The two data cards contain alphanumeric identifiers including the names or codes of the 
aircraft whose P s data are on the tapes. The program operates as follows. The data 
tapes created by Program AMP are searched for the P s value corresponding to a given 
altitude, Mach number, and normal load factor. Next, the turn rate and turn radius for 
the specified Mach number and normal load factor are computed. The Pq values found 
on the data tape and the corresponding turn radii and turn rates are then used as coordi- 
nates for one of the Langley graphics routines. 

In the current version of NEMO, nine plots are developed. These plots are the 
resulting combinations of three altitude levels and three Mach number levels. The alti- 
tude levels are 10 000 feet (3.05 kilometers), 20 000 feet (6.1 kilometers), and 30 000 feet 
(9.14 kilometers). The Mach levels are 0.6, 0.9, and 1.1. The foregoing values of altitude 
and Mach number are those used in the NASA/Navy study but they are variable and any 
Mach-altitude set can be developed. A typical plot resulting from NEMO is given in fig- 
ure 5. As shown in the figure, the variations of P s with turn rate for both aircraft are 
plotted on the same set of axes for comparative purposes. Symbols along these curves 
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correspond to integer levels of normal load factor as indicated. Also plotted in this fig- 
ure is the corresponding turn radius associated with the turn rate. 

Program EXTRA: Program EXTRA takes the P s data for two aircraft (furnished 
by two AMP tapes) and computes the difference in specific excess power AP S between 
the two aircraft at various Mach numbers, altitudes, and normal load factor levels. 

These AP S values are then used in various forms of output. 

The input to the program consists of two AMP tapes and aircraft identification 
cards. The data are taken from the tape and placed in arrays as a function of Mach num- 
ber and altitude. Values of AP S are then computed by using these arrays. At this 
point, the AP g data are used in an iteration process similar to that used in Program 
CONTOUR. This consists of varying Mach number and altitude for a given load factor 
to determine Mach number-altitude sets for a selected value of AP S . This gives infor- 
mation necessary for AP S contour plots. 

Two types of output are provided with the AP S contour information. One is a plot 
tape which can be used in one of the Langley plotting systems. This option has not been 
used. The other output consists of line printer listings of the AP S data for selected 
Mach numbers, load factors, and altitudes, and a line printer plot of AP S superiority 
regions for the aircraft combination being considered. Figure 6 illustrates the line 
printer type of plot. As shown, the plot lists the altitude and Mach number conditions 
where either aircraft has a AP S advantage or superiority of 50, 100, or 150 feet per 
second (15.2, 30.5, or 45.7 meters per second). Separate plots are generated for a num- 
ber of load factors. The AP S map represented by figure 6 can be used to show areas 
where one aircraft has a potential advantage in energy and thus where the pilot should try 
to maneuver it during the ACM. 

Program SUGMUG: Program SUGMUG takes the same DMS aerodynamic input 
data (for a single aircraft) that were used in Program AMP. The program determines 
numerical values for sustained and maximum load factors for various combinations of 
Mach number and altitude. 

The program uses three loops: altitude, Mach number, and angle of attack. It 
operates by selecting an altitude and a Mach number and iterates angle of attack through 
a range of values. At each value of a, a value of lift coefficient, drag coefficient, and 
thrust is obtained. A check is performed to ascertain that the value of C L does not 
exceed Cl max- Computations are performed to determine n z and P s for this set 
of conditions. Angle of attack is iterated until P s becomes zero or negative or maxi- 
mum lift is exceeded. If it is zero, then the values of Mach number, angle of attack, nor- 
mal load factor, and altitude are placed in a holding array for later use. If P s is not 
zero, but is negative, then it is used with the previous positive value of P g to interpo- 
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late and obtain values of angle of attack, normal load factor, Mach number, and altitude 
at P g = 0 and the values are placed in the holding arrays. If P s is positive, the iter- 
ation of angle of attack continues. With this technique, values of a and n z are obtained 
when P g = 0 for a specified set of altitudes and Mach numbers. 

This same type of process is employed for the maximum load factor case. However, 
instead of doing all the calculations to obtain P s = 0, the calculations are performed by 
using the maximum lift coefficient obtainable for each Mach number and altitude combi- 
nation. From these calculations, values of n z , a, and P s are obtained at maximum 
load factor for a specified set of Mach numbers at a number of different altitudes. 

Data for both sustained and maximum load factors are generated as line printer out- 
put and are also placed on a magnetic tape for use in Program UNICORN. 


Program UNICORN: Program UNICORN takes up to three SUGMUG tapes and pre- 
pares load-factor and angle -of -attack data from the tapes for plotting. 


The input to UNICORN is the SUGMUG tape, or tapes, and labeling information for 
resultant plots. The program takes the SUGMUG data and determines the locations on 
a plot and stores the vertical and horizontal locations on a plot tape. This operation is 
done for both n z and a? as a function of Mach number. The resultant plot tape is used 
in one of the Langley plotting systems to give plots of n z and a as a function of M 
for each altitude and for sustained and maximum n z . A typical output of UNICORN is 
given in figure 7 for two aircraft. This plot represents a sustained n z case. 


Program SVELTE: Program SVELTE creates a plot file of maximum specific 
excess power P g max as a function of specific energy E s at given levels of turn 
rate for a number of aircraft. When more than one aircraft is being considered, plots 
of differential maximum specific excess power AP g max are also developed. 


The input to the program consists of the same DMS performance data packages 
used in AMP and SUGMUG and cards with alphanumeric identification. The program 
operates as follows. Values of turn rate and E s are specified. Mach number starts 
at zero and is increased as altitude is decreased with E Q constant 


E s = h + If) 


At 


each turn rate, specific energy, and altitude, the velocity, Mach number, normal load fac- 
tor, and P g are calculated and the P S)ma x over the altitude range is updated. When 
the altitude decrementation is completed, the E g value is placed in the abscissa array 
and the P S; max value is placed in the ordinate array, and these data are placed on a 
plot file. A new E s value is then specified and the ordinate calculation is repeated. 
This process is continued until the desired number of turn rate values are examined. 


The plot file is the primary output of Program SVELTE. This plot file is used to 
provide curves of P s max as a function of E g . An example of SVELTE output is given 
in figure 8. 


10 



A P s max plot consists of a number of curves of P S}max as a function of E s , 
each curve corresponding to a different turn rate for one aircraft. After a P s max plot 
has been developed, plots of AP s rnax between the present aircraft and each one of the 
previously considered aircraft are developed. The entire process is repeated for each 
aircraft. In addition to the plot file, similar data are output on the line printer. This 
information includes the abscissa and ordinate (E g and P S)max or AP s max j of each 
point plotted, the value of turn rate at which each point is plotted, and the, altitude, Mach 
number and normal load factor at which maximum specific excess power occurred. 

Daily Use Programs 

To meet requirements for daily single run data reduction, display, and analysis, 
one primary and several secondary data reduction and analysis programs have been gen- 
erated. As mentioned previously, the daily output of these programs was used by pilots 
and test engineers in the NASA/NAVY study to evaluate the previous day's results. The 
results were evaluated to determine whether maneuvers used were successful or not and 
whether any unusual outcome was encountered. In addition, these data provided the basis 
for the final evaluation of the tests. The programs that were used and that are available 
are shown in figure 9. As can be seen in the figure, the primary input to all the programs 
is the DMS Data Tape which was discussed previously. Each of these programs and their 
elements are discussed herein. 

A general comment related to all programs shown in figure 9 is the heading infor- 
mation present on all output. During the DMS runs, information was placed on the DMS 
Data Tape in the integer word KODE to facilitate record keeping and future use. KODE 
contains the date, run number, aircraft pair being flown, pilots in each aircraft, and the 
I.C. for a given run. In each of the following programs, KODE was used in identifying 
runs and as information for output headings . 

Program HASTE .- Program HASTE was developed by Naval Ship Research and 
Development Center (NSRDC) personnel to perform the major data reduction and display 
functions required to support the NASA/Navy study. The primary operations of HASTE 
include reading the DMS Data Tape, computing variables of interest from the DMS vari- 
ables, and calling and controlling data handling subroutines. Figure 10 shows the data 
handling subroutines that can be called by HASTE and the type of output they produce. As 
seen in the figure, the subroutines and the output they produce are as follows: 

(1) Time history listings of 68 variables by subroutine TLIST 

(2) Time history listing of 30 event markers by FLAGER 

(3) Cumulative histograms of 10 variables from FASTON 

(4) 24 two- and three-way histograms by TTHIS 
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(5) Time history graphical display of 21 variables for up to 900 seconds of flight 

time and two X,Y-type plots every 300 seconds by GRAFIC 

(6) Gun fire and offensive position summary section by S1RUN 

(7) The first, second, and third cumulative central moments of 10 variables at 

various times during a DMS run presented by M123S 

The output indicated are discussed in the descriptions of the subroutines. 

The input to HASTE consists of the DMS Data Tape and a HASTE Data Package. 

The HASTE Data Package includes the following information: 

(1) Gun and aircraft vulnerability definition 

(a) Number and types of guns 

(b) Characteristics of each gun type 

(c) Vulnerability profile (target area at different ranges, altitudes, etc.) 

(2) Number of runs to be processed 

(3) Entry rates for subroutines (How often is each subroutine to be entered ?) 

(4) What subroutines are to be used (What type output is desired?) 

(5) Frequency of histogram, data listing, etc., output 

Program HASTE operates in the following manner. First, the HASTE Data Package 
is read and all variables used in HASTE are initialized. Next, each DMS Data Tape file 
(run) processed by HASTE is read one record at a time. This is followed by computation 
of the variables required for the subroutines. Then, a test of the indicators associated 
with each subroutine is executed and the necessary subroutines are called. After all the 
subroutine calls, another DMS record is read and the computation, testing, and calling 
cycle is repeated. This cycle is continued until a run is processed. When the run has 
been completed, every subroutine that has been used is called again for final disposition 
of its accumulated data. This process is continued until all runs have been completed. 

The output resulting from HASTE is dependent on which of the optional data handling sub- 
routines shown in figure 10 are called for in the HASTE Data Package. The data subrou- 
tines are independent of each other and under programer control in that data in the input 
card deck prescribe the frequency of display of each display and the time between samples 
with which the displays receive data. The fact that these values are supplied by the pro- 
gramer provides, first of all, flexibility to process data tapes which have radically differ- 
ent sampling rates. Secondly, programer control over the displays provides the ability 
to tailor-make the output to analysis needs with regard to type, amounts, and accuracy of 
displays. A brief description of the subroutines and a typical output of each follows. 
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Subroutine TLIST: Subroutine TLIST accumulates data from the DMS Data Tape 
and HASTE and outputs time listings of 68 parameters. Three types of output are pro- 
vided, with one page per aircraft for each type (six pages per specified interval of time). 
The first two pages contain time histories of aircraft relative positions and velocities. 

The parameters given include time into run t, range between aircraft R, the time rate 
of change of range R, the difference between the inertial coordinates of the two aircraft 
AX, AY, AZ, the difference between the inertial velocity components of the two aircraft 
AX, AY, 4 AZ, and the angular rate of change of the range vector in the inertial coordi- 
nate system d,p and its components in the vertical (X-Z plane) and horizontal (X-Y plane) 
directions dy, djj. The next two pages contain aircraft specific information. The 
information presented consists of each aircraft's altitude h; velocity V; Mach number 
M; Euler roll, pitch, and yaw angles cp, 9, ip and angular velocities p, q, r; normal 
and lateral load factors n z , ny; angle of attack a\ sideslip angle j3; specific excess 
power P g ; and an AML cell score similar to that generated in the AML program previ- 
ously discussed. The last two pages of the series contain time listings of angular rela- 
tions between the two aircraft. This type of output includes the angular relations illus- 
trated in figure 11. The relations indicated in this figure are the off-boresight angle X 
which is the angle between the x-axis and the range vector R, the deviation angle e, 
which is the angle between the velocity vector V and range vector R. The other angle 
given is angle-off. Angle-off of an attacker aircraft is defined as the difference between 
the target aircraft's deviation angle and 180°. In addition to these angles, various com- 
ponents of X are also presented (X u y D , X w , X e , xA 

Each type of information presented by TLIST is output for both aircraft at user 
specified time intervals. Examples of the three types of TLIST output are given in 
tables V to VH. 

Subroutine FLAGER: Subroutine FLAGER takes the flag word (IFLAG) from the 
DMS Data Tape, decodes it, and uses values contained in it to generate a page of event 
markers as illustrated by table VIII. The flag word is an 8-digit integer, with each digit 
representing one event and its status. The events contained in IFLAG and shown in 
table Vin are as follows: 

(1) Gun status - Is aircraft in gun envelope (E) and/or firing guns (F)? 

(2) Missile A status - Is aircraft in missile A launch envelope (simplified launch 
cone defined by X and R) (E) ? Is missile A selected (S) ? Is missile A launched 
(F)? 

(3) Repeat for missile B. 

(4) Repeat for missile C. 
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(5) Bingo marker - Has aircraft reached bingo fuel limit or not (bingo fuel is min- 
imum fuel required to return to safe area) ? 

(6) Angle-of-attack marker — Indicates when angle of attack exceeds designated 
maximum value. 

(7) High lift device activation marker - Indicates when high lift devices are used 
(flaps, slats, etc.). 

(8) Ground impact marker - Indicates when aircraft reaches altitude less than zero. 

These events are given for both aircraft. The run shown in table VIH indicates, for exam- 
ple, that sphere A had missile type A selected and was in the envelopes of both missile 
types A and B from 160 to 170 seconds. At 170 seconds, sphere A fired a missile type A 
and then went out of the missile envelope for the rest of this page. Sphere A also had an 
angle of attack exceeding 30° at 160 to 162 seconds and later at 186 to 188 seconds into 
the run. Sphere B had missile B selected from 160 seconds until 180 seconds. During 
that time sphere B fired a missile type B, while in envelope, at 178 seconds. At 184 sec- 
onds sphere B selected guns and fired them, while in envelope, for the period between 
184 and 196 seconds. The time interval for printing the FLAGER page is an input 
parameter. 

Subroutine TTHIS: Subroutine TTHIS accumulates data from the DMS Data Tape 
and from HASTE and uses them to compute and write values for time cumulative two- 
and three-way histograms. Table IX lists the 24 histograms provided by TTHIS. A 
representative three-way histogram is shown in table X. This table gives the percent- 
age of the total run time within various interval sets of the off-boresight angle X for 
aircraft A and B and the range between the two aircraft. Table X represents only part 
of this particular histogram since X for aircraft A is shown only up to 60°. This 
particular histogram can be used to determine the percent of time one aircraft spends 
in the envelopes of various weapons. The histograms provided by TTHIS can be output 
at any interval desired by the user through HASTE Data Package input. An example of 
this option would be an output every 60 seconds up to the total run time, that is, at 60, 

120, 180, and 240 seconds. 

Subroutine FASTON: Subroutine FASTON computes and displays the percentage of 
time in which selected parameters, indicated in table XI, lie within chosen intervals. 

These percentages are cumulative with respect to time into flight and computed at equally 
spaced time values. An example of one of the displays is given in table XH. This table 
gives the percent of time spent in various intervals of altitude as a function of time into 
the run (every 25 seconds up to 225 seconds). These data are presented for both aircraft. 
The time into run intervals are varied through the input to the HASTE Data Package. 
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Subroutine GRAFIC: Subroutine GRAFIC takes 21 parameters from HASTE and 
prepares them for graphical display by one of the Langley graphics systems. When the 
data for a run are accumulated, these accumulated data are placed on a magnetic tape 
which is used to provide graphic output. The output provided by GRAFIC consists of 
nine frames of graphics per run. 

The nine frames of an output set consist of two different types of plots; those which 
plot a parameter value against time and those which plot one parameter's value at a time 
against another's value at that same time. 

The first seven frames plot 21 parameters against time. Frame 1 has two plots 
and is illustrated in figure 12. On the lower plot, range is plotted as a function of time; 
on the upper plot, off-boresight angle is plotted as a function of time for both aircraft A 
and B. Off-boresight angle for aircraft A is plotted with + signs and off-boresight angle 
for aircraft B, with dots. The other parameters, all plotted against time, are as follows: 


Frame 2 

R, 6 


Frame 3 

h(l), V(l), h(2), 

V(2) 

Frame 4 

n z (l), a( 1), 0(1), 

P S (1) 

Frame 5 

n z (2), Q!(2), 0(2), 

P S (2) 

Frame 6 

e(l), Angle- off (1) 


Frame 7 

e(2), Angle -off (2). 



The frequency with which points are plotted along the time axis is controlled by user 
input data. 

The last two frames of data output by GRAFIC plot altitude against Mach number. 
The plot gives Mach-altitude points for every 10 seconds into the run and indicates by 
symbol whether at that time the aircraft was on offense (A), defense (X), or in a neu- 
tral position (N) with respect to other aircraft. This output is illustrated by figure 13. 
A frame is provided for each aircraft. An aircraft was defined as on offense when its 
A was less than 90° and the other aircraft's A was greater than 90°. An aircraft was 
considered to be on defense when the other aircraft was on offense. A neutral position 
was indicated when neither aircraft was on offense. 

Subroutine S1RUN: Subroutine S1RUN compiles statistics and histograms which 
summarize gun and missile zone and firing information and offensive position data. The 
output from SlRUN is presented on six different pages. One of the prime relations con- 
sidered is the off-boresight angles A and the range between the two aircraft. This set 
of information indicates the offensive status of each aircraft and whether they are within 
simplified weapon envelopes. The pages output and what they present are as follows. 
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The first page of output is illustrated by table Xin and includes the following 
information: 

(1) Time of first entry into a given weapon envelope (guns and three missiles) 

(2) Total time on offense (X of attacker < 90°) as a function of time into run 

(3) Total time on offense with advantage (x of attacker < 90° and X of tar- 

get > 90°) as a function of time into run 

(4) Total time in gunfire envelope (100 feet < range < 3000 feet, target X > 120° 

and attacker X < 1°, 2°, 3°, 5°, and 10°) as a function of time into run 

(5) Same as (4) except that target X > 90° 

(6) Total time guns were firing while in envelope for envelopes given in (4) and (5) 

(7) Time in each of the three missile envelopes as a function of time 

The second page of output is illustrated by table XIV and gives percent of gunfire time 
within one-way intervals of deviation angle, range, and angle -off. The third and fourth 
pages of S1RUN output are similar and are represented by table XV. These pages give 
percentage of offensive time and offensive time with advantage, respectively, within one- 
way intervals of range, range rate, altitude, deviation angle, Mach number, and normal 
load factor. The fifth page of output is illustrated by table XVI and gives a listing of each 
gunfire burst and missile launch initiated by a DMS pilot during a run. The gunfire output 
includes the start and stop times of the burst and the probability of kill (APK) of the gun- 
fire burst. This page also has the total run time, the accumulative gun APK, and the 
average AML score for a run. The missile launch output, similar to gunfire output, gives 
an indication of each DMS trigger pull when a missile had been selected, the type of mis- 
sile selected, the time of the trigger pull, whether or not the missile was launched in 
envelope, and the aircraft the missile was launched from. There is no evaluation of mis- 
sile probability of kill P^ in HASTE, and no firing logic, weapons load considerations, 
etc., are indicated in these data. 

The five pages of output given by S1RUN were used as the primary method of daily 
run review by the test pilots and were generated for all data runs. In addition to the line 
printer output from SI RUN, the first four pages are output on punched cards or magnetic 
tape for use in Program MULTI2. 

Subroutine M123S: Subroutine M123S takes selected DMS data (R, R, V, h, n z , 
X) and calculates and lists the cumulative first and second central moments at given times 
into the run. The first central moment corresponds to the means of the parameters over 
the specified time interval and the second moment corresponds to the standard deviations. 
Subroutine M123S also outputs the initial conditions of X a , V, h, cp, 9, \p, and R 
for each run. Typical M123S output is given in tables XVII and XVM. 
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NWC missile program .- The NWC missile program is one of two off-line mis- 
sile programs used to evaluate weapons, weapon loads, and launch opportunities and suc- 
cess. It was developed by personnel of the Naval Weapons Center (NWC) and modified 
by NSRDC and LRC personnel for the DMS. The NWC missile program utilizes aerody- 
namic descriptions of current missiles. These descriptions are used by the program to 
compute missile trajectories after launch. The trajectories are continued until impact 
or miss, due to exceeding some missile limit, occurs. Various output options are 
available . 

The input to the program includes the DMS Data Tape and a data package which _ 
describes missile aerodynamic characteristics, thrust profiles, sensor gimbal angles 
and rates, and various limits (launch envelopes, sensor tracking rates, miss distance 
criteria, etc.). Any number of missile types can be evaluated in the program although 
only two types were used in the NASA/Navy study. 

The program operates by investigating the situation of one aircraft at a time, with 
only one type of missile being evaluated at a time. The program takes the DMS Data 
Tape and checks at each time interval on the tape whether the aircraft meets the launch 
criteria and if not what was the reason. The launch criteria consist of off-boresight 
angles, ranges, line-of-sight sweep rates, etc. If the criteria were met, a launch is ini- 
tiated and the resulting trajectory is calculated. If an impact or a miss within a specified 
distance occurs, the program indicates an impact and various conditions at impact. If 
the missile exceeds a limit during its flight, the trajectory calculations are discontinued 
and the reason is indicated. If an aircraft does not satisfy the launch envelope criteria, 
the reason is indicated and no trajectory is calculated. This check occurs for each time 
interval on the tape or at an interval specified by user. The launch criteria checks and 
successful launches are counted during a run. All runs on a tape are evaluated for the 
particular missile and aircraft combination and then the DMS Data Tape is rewound and 
the run evaluated from the other aircraft's viewpoint. This process continues until all 
aircraft and missile combinations have been evaluated. 

The output of the NWC missile program is illustrated in tables XIX and XX. 

Table XIX is one page of time listings for a run. Below the standard heading is a line 
of print which gives the code words for the missile and aircraft combination being eval- 
uated. Next is a line indicating the launch criteria for the missile being evaluated. These 
criteria include the launch gimbal angle limit (missile sensor gimbal angle), the maximum 
and minimum range limits, the sensor track rate limit, and the maximum gimbal angle 
limit after launch. The rest of the data given in the table relate to launch attempts. The 
first section, labeled "Launch Conditions," gives the conditions at selected time intervals 
when a launch is attempted. This section includes the range, sensor gimbal angle, and the 
sensor tracking rates. If any of these values exceed the launch limits, then a launch can- 
not occur and a statement telling which envelope criteria were not met is printed. These 
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statements are illustrated by the printouts for the time interval between 108.5 and 
111.0 seconds. Next is given the values which indicate a stop after a launch. This sec- 
tion is labeled "Limit Stop." This section includes the time interval after launch a limit 
stop is exceeded (T2), the gimbal angle (GA2), the sensor tracking rate (TR2), and the 
stop limit range (SLRNG). If any of these conditions exceed the limit criteria, then the 
trajectory is stopped and the reason printed out. These reasons are illustrated by the 
printouts for the period between 104.0 and 108 seconds. Finally, if a launch continues on 
to impact, the final section, labeled "Intercept Data," is printed out. This section gives 
conditions at impact or closest approach and maximum values occurring during trajectory 
The values given include the maximum load factor (GN) and angle of attack (ALM) 
occurring during trajectory, time of missile flight (TOF), miss distance (MSD), clos- 
ing velocity (CLOVEL), range at impact (RSTP), yaw, pitch, and roll angles (AANG, 
EANG, RANG), missile velocity (MVEL), cross angle (CANG), angle of attack (ALF), 
off-boresight angle (LANG), longitudinal acceleration (GAX) at impact, and time of 
impact (TIMP). This output is illustrated by the line printed at 111.5 seconds in 
table XIX. In addition to the programed attempts to launch a missile, any trigger pulls 
recorded during a DMS run, with the correct missile selected, will cause a special print- 
out indicating the launch and its success or failure. Table XX is an example of the 
summary-type page presented after each run. This output tabulates the simulated and 
DMS launches and their successes. The first set of data indicates the number of simu- 
lated launches out of envelope, in envelope but missing because of a limit stop, and hits. 
The frequency of launch attempts is also indicated. A similar output is given for pilot 
initiated launches. The final set of data incorporates a time in envelope delay and counts 
only the launches meeting this additional constraint. This output presents the total num- 
ber of launches with the delay and the number of impacts with gimbal angles at launch of 
less than or equal to 40° and 20°. 

i 

NADC missile program (LAUNCH) .- The NADC missile program was developed for j 
the Naval Air Development Center (NADC) to support the NASA/Navy DMS study. This 
missile program utilized launch acceptance region (LAR) tables, based on actual results 
in combat and in extensive simulations, to determine when a successful launch could be 
initiated. The LAR tables contain combinations of range, off-boresight angle, target 
aspect angle (angle between target x velocity axis and range vector), Mach number, and 
altitude which should provide a successful missile impact. A graphical representation 
of a part of the LAR table for one missile is given in figure 14. As shown, this table is 
for the Mach 0.9 region, an altitude of 15 000 feet (4.6 kilometers), heading angles of 0° 
and 30°, and an assumed target maximum acceleration of 6.5g. When conditions within 
these acceptable zones are met, a successful launch is assumed. Various inputs are 
provided to describe the results of a run. 
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Inputs to the program consist of the DMS Data Tape and the LAR table information 
for up to three weapon systems in the form of data cards. These cards include missile 
launch sequencing and logic and run information indicating runs to be processed and wea- 
pon systems on each aircraft. 

The program operates in the following manner: 

(1) Reads LAR table and weapon constraint input cards 

(2) Reads run information cards 

(3) Reads DMS Data Tape record for current run 

(4) Initializes run variables and arrays 

(5) Computes states of both aircraft for current record 

(6) Computes launch zone conditions for all available weapons 

(7) Tests launch criteria for all available weapons 

(8) Determines impact conditions for all launched missiles 

(9) Prints launch report data for current record 

(10) Continues (1) to (9) until run completed 

(11) Outputs run summary data 

The output provided by LAUNCH is illustrated in tables XXI and XXII. Table XXI 
is one page of the launch. report output for each time interval on the DMS Data Tape. Each 
page of the launch report has a heading describing all the run information including the 
code for the weapons on each aircraft. Below this heading, the launch data for both air- 
craft are listed. For each time interval each aircraft’s altitude, Mach number, normal 
load factor, off-boresight angle, aspect angle, and range are given. In addition, for the 
missiles carried, indications of when the aircraft are in a missile launch zone are listed. 
The letter H indicates that an aircraft has a proper heading but is not in the launch zone 
with respect to aspect angle or range. The letter V indicates that an aircraft is in the 
launch zone with respect to all parameters. The other columns indicate simulated and 
pilot controlled launches with the letter V representing simulated launches when time 
in zone criteria are satisfied and dashes indicating impacts for both the simulated and 
piloted launches. In addition, if a pilot attempts a launch while out of the envelope, the 
letter T is used. 

Table XXII shows the launch summary report and impact matrix output at the end 
of a run. The launch summary report gives the DMS Data Tape numbers, run numbers, 
date, pilot, aircraft codes, weapons codes, number of time intervals in launch zone, and 
simulated launches for each missile type for each aircraft. The impact matrix outputs 
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the simulated and piloted launches, indicating the type launched and the time at which 
it occurs. The code for the type of missiles launched is also given. In addition to the 
line printer output, the impact matrix and launch summary report are output on punched 
cards for use in Program OUTCOME which is discussed in the section entitled "Posttest 
Programs." 

Program 3-D .- Program 3-D takes data from the DMS Data Tape and generates 
plots which represent a ground trace (X-Y plane) and a perspective 3-D representation 
of the two aircraft and their trajectories as seen from a selected point in space. The 
only input to the program is the DMS Data Tape and some identification and scale factor 
cards. The program operates in the following manner: 

(1) The inertial coordinates of the center of gravity and the attitude of each aircraft 
are read from the Data Tape. 

(2) Each aircraft is assumed to be represented by a simple nine-point vector figure, 
defined in the body axes. The computer uses the aircraft attitudes to transform the nine 
points from the body axis system to the inertial coordinate system. 

(3) By using a predefined eye position and point of regard, the coordinates of the 
points are transformed from the inertial system to a line-of- sight coordinate system 
having the X-axis alined with the line of sight. 

(4) By assuming that the plane of the plot is scaled to a location at some distance 
along and perpendicular to the line of sight, each of the points is located on the 3-D plot. 

(5) The nine points are connected to form the aircraft symbol 

(6) The locations of the two points representing the wing tips are stored for use in 
succeeding plots to define the trajectory of the aircraft. 

(7) The velocity, range, and time into the run of each aircraft are printed. 

The only difference hetween the 3-D perspective plot and the ground trace plot is 
that it is assumed for the ground trace that the aircraft is at zero altitude (Z =0). Fig- 
ure 15 is representative of the output of Program 3-D. In the figure the solid lines repre- 
sent one aircraft's trajectory and the dotted lines represent the other aircraft's trajec- 
tory. One pair of side-by-side dots represents the aircraft wing tips and these data are 
output at a specified time interval. In this figure the data are output every second with 
20 seconds given on each frame. The aircraft symbols and alphanumerics are given 
for the last time on the frame. 


Posttest Programs 

After a series of runs were completed for an aircraft pair, it was a requirement of 
the NASA/Navy study to perform multirun analysis and data presentation. This type of 
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information gave some statistical support to trends encountered in single runs and per- 
mitted more concise presentation of study results. Three programs were used as repre- 
sented by the dashed boxes in figure 9. All the programs are generally concerned with 
the aircraft's weapons and offensive performance with two of them being concerned pri- 
marily with time in zones and the other one being concerned with weapons performance 
and probability of kill P^. Two of the programs output single run data as well as cumu- 
lative or average values. A discussion of these three programs follows. 

Program OUTCOME .- Program OUTCOME was developed for NADC and its pur- 
pose is to determine the probabilities of win, lose, or draw in a simulated combat between 
two aircraft. Program OUTCOME receives its primary input data from the punched out- 
put of the LAUNCH program. Program OUTCOME determines sets of win, lose, or draw 
probabilities based upon sequences of both simulated (i.e., possible) impacts and actual 
(pilot initiated) impacts coming from LAUNCH and missile probability of kill and weapon 
load limitation parameters input by the user. In general, the probabilities of outcome 
are reported in a matrix form which is the result of a parametric study. 

Program OUTCOME accepts the probabilities of kill for the firing of a single mis- 
sile from each weapon system. It then varies the probability of kill of one weapon sys- 
tem (e.g., missile type) in steps from 0 to 1 and observes the effect upon the probabili- 
ties of win, lose, or draw. The program next varies a second probability of kill (either 
another weapon system on the same aircraft or a weapon system on the other aircraft) 
against the first while holding the remaining probabilities of kill constant. The variation 
of one probability of kill with another creates the matrices of win, lose, or draw. 

For each run the program computes and prints matrices of win, lose, or draw for 
both the "simulated" and the "piloted" cases. At the end of the program, matrices which 
are the averages of the individual run matrices are computed and printed. An example of 
the output from OUTCOME is given in table XXIH. This table presents an average prob- 
ability of loss matrix for sphere A for a set of runs. The first set of information on the 
left side of the page indicates what weapons are to be considered on each aircraft and the 
number of each type. Below that is given the desired single shot probability of kill for ' 
each weapon. The missiles with = 1.0 are those to be varied from 0 to 1.0 in the 
matrix while the other missile values are to be held constant during computation. 
Other heading information given by this output includes the DMS Data Tape being pro- 
cessed, the run number, the time interval in the run to be evaluated, and the code for 
the weapons that are being evaluated in the probability matrix. The matrix size itself 
can be varied by the user. 

Program MUL-TI2 .- Program MULTI2 takes run data output by subroutine S1RUN 
of Program HASTE and combines it according to pilot, date, I.C., or aircraft as specified 
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by input cards. These data are then used to determine means, standard deviations, and 
t-test results for a specified set of runs and two histograms. 

Input to MULTI2 consists of the punched cards or magnetic tape representing single 
run results coming from HASTE, cards which identify parameters that will be considered 
and output, and cards that identify runs to be used in computations. One set of run cards, 
including an alphanumeric identifier, is necessary for each set of data to be analyzed. 

Runs to be evaluated as a set are designated by user. For example, if it is desired 
to evaluate the performance of one pilot against all other pilots, the runs in which the 
selected pilot was involved is specified as a set according to run dates. Any other eval- 
uation would require another set of cards with desired runs listed. 

The program operates in the following manner. First, it reads the alpha-numeric 
title for the set of data being processed, determines the number of parameters to be eval- 
uated, and uses the number of parameters to establish a set of holding arrays. Next, 
MULTI2 reads and stores in arrays dates and runs to be used in the evaluation. When 
the run and date arrays are established, the program reads the run data from HASTE 
and checks each run number on tape with those in array. When a date matches, the 
information on that run is stored in the proper parameter arrays. This process con- 
tinues until all runs have been located. The parameter data are then sent to output sub- 
routines which perform the various computations required to obtain the means, standard 
deviations, t-test parameters, and histogram information. These data are printed out, 
the arrays are initialized, the next set of runs is read in, and the process is repeated 
until all sets of runs are processed. 

Typical MULTI2 output is shown in tables XXIV to XXVIII. Table XXIV is the same 
output as given in table XM except that the numbers shown represent mean values for a 
set of runs instead of one run. The heading describes the particular set of runs being 
considered. Table XXV presents also the same parameters as discussed in table XIII 
except that the numbers represent standard deviations. Table XXVI contains the values 
of the t-test (ref. 5) and the number of degrees of freedom (D.F.) for use in testing the 
hypothesis that the mean values for the two aircraft (table XXIV) are equal. Table XXVII 
gives a three-way histogram similar to one output by TTHIS of Program HASTE which 
gives the mean percentages of time at various intervals of n z , M, and h. Table XXVin 
is a three-way histogram similar to that discussed in table X except that it too represents 
mean values for a set of runs. 

Program BZONE .- Program BZONE was developed to evaluate a pilot's capability 
to establish and maintain position of one aircraft with respect to another. This evaluation 
is accomplished by establishing three target zones and five attacker off-boresight-angle 
zones. The three target zones are illustrated in figure 16. The three target zones cor- 
respond to a head-on, beam, or tail situation. The five attacker zones correspond to off- 
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boresight angles of ilO°, §20°, ^40°, ^60°, and <90°. An additional restriction of range 
of <3000 feet (0.914 kilometer) is included for the =10° zone. The resulting 15 attacker- 
target zone combinations give general trends as to relative offensive superiority and 
weapons delivery opportunities. 

The only input to BZONE is the DMS Data Tape. Program BZONE operates with 
these data by reading the inertial positions and Euler angles and determining both air- 
craft's off-boresight angles. The angles, in turn, are checked to determine into which 
one of the 15 combined zones they fall. When the zone for the current record has been 
determined, the program checks to see whether it is the first time the zone has been 
entered and updates a time in zone counter. This process continues until a run is com- 
pleted. After a single run is processed, the results are output and the first entry times 
and total time in zones are transferred to a set of holding arrays which are used for mul- 
tirun averaging. The rest of the runs are processed in a similar manner until all runs 
are completed. At this time the holding arrays are used to compute average values of 
all the single run parameters. A single page of average values are output for each I.C. 

Typical output from BZONE is given in tables XXIX and XXX. Table XXIX is repre- 
sentative of a single run output. The headings across the page relate to the target air- 
craft's three orientations as shown in figure 16, and the vertical headings are the 
attacker's off-boresight-angle intervals. The first three columns show the times of 
first entry into each of the 15 zones. The next three columns show total time in each 
of the zones. The last three sets of three columns indicate the first time in which the 
attacker aircraft entered one of the zones and remained in it for 2, 4, and 6 seconds, 
respectively. Directly below these columns are corresponding columns which indicate 
the number of times during the run in which the 2, 4, or 6 seconds in zone occurs. The 
same type of data is also given with the attacker and target roles reversed. 

Table XXX represents the multirun output for BZONE. The same type of headings 
are used for the multirun output with the numbers given representing average values 
instead of single run values. The averages given, however, are not determined by the 
total number of runs being considered but instead are based on the total number of runs 
in which a zone has a nonzero value; that is, if one run had no first entries into one of 
the zones, that run would not be counted in averaging the first entry times for a set of 
data. The number of runs with nonzero values for each zone and for each element of 
output is output directly below the element of interest. The number of runs also indi- 
cates the degree of reliance that can be placed on average values, when compared with 
total number of runs. The output is again given with both aircraft as attacker and 
defender. 
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CONCLUDING REMARKS 


An extensive data reduction and analysis capability has been developed to support 
differential maneuvering simulator (DMS) aircraft simulation studies. This capability 
includes programs which provide (1) pretest information which can be used to show 
pilots' areas of aircraft superiority and possible maneuvers, (2) daily use information 
which presents individual run data in different forms, and (3) multirun data to show, con- 
cisely, results of a series of runs. 

These programs use, in general, one of three types of input data depending on which 
one of the three types of program is involved. The pretest programs utilize as input the 
same aerodynamic data used in the DMS to describe the aircraft of interest. Daily use 
programs, as described, use output from the DMS simulation as the major input. Multi- 
run programs, in general, use as input the output from the daily use programs. The 
major programs are designed to be flexible and are generally adaptable to any aircraft 
pair . 

The discussion of each program included in this paper was not intended to be a 
user's guide but to be an introduction to each available program and major input and 
output data. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 5, 1973. 
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TABLE I.- DATA STORED ON DIFFERENTIAL MANEUVERING SIMULATOR DATA TAPE 


Word 

Variable 

Description 

Unit 

1 

KODE 

Code number of date, run, pilot, and aircraft 

- - - 

2 

T 

Time from start of run 

Seconds 

3 

RBA (R) 

Range 

Feet 

4 

RBADOT (R) 

Range rate 

Feet/second 

5 

IFLAG(l) 

Integer FLAG for gun/missile data for sphere A 

- - - 

6 

IFLAG(2) 

Integer FLAG for gun/missile data for sphere B 

- - - 

7 

X(l) 

X inertial position of sphere A 

Feet 

8 

Y(l) 

Y inertial position of sphere A 

Feet 

9 

HAB(l) (h) 

Altitude of sphere A 

Feet 

10 

SXDOT(l) (X) 

X inertial velocity of sphere A 

Feet/second 

11 

SYDOT(l) (Y) 

Y inertial velocity of sphere A 

Feet/second 

12 

HDOT(l) (h) 

Time rate of change of altitude of sphere A 

Feet/second 

13 

PSID(l) W 

Euler yaw angle of sphere A 

Degrees 

14 

THETAD(l) ( 9 ) 

Euler pitch angle of sphere A 

Degrees 

15 

PHID(l) ( <p ) 

Euler roll angle of sphere A 

Degrees 

16 

NZ(1) 

Normal load factor of sphere A 

G’s 

17 

NY(1) 

Lateral load factor of sphere A 

G's 

18 

ALPDEG(l) (a) 

Angle of attack of sphere A 

Degrees 

19 

BETADEG(l) (/3) 

Angle of sideslip of sphere A 

Degrees 

20 

P(l) 

Body roll rate of sphere A 

Radians/second 

21 

Q(D 

Body pitch rate of sphere A 

Radians/second 

22 

R(l) 

Body yaw rate of sphere A 

Radians/second 

23 

WEIGHT(l) 

Weight of sphere A 

Pounds 

24 

TT(1) 

Thrust of sphere A 

Pounds 

25 

PS(1) 

Specific excess power of sphere A 

Feet/second 

26 

DELTSP(l) 

Fore/aft stick position of sphere A 

- - - 

27 

DELTAP(l) 

Right/left stick position of sphere A 

- - - 

28 

DELTRP(l) . 

Rudder pedal position of sphere A 

- - - 

29 

TPOSTOL(l) 

Throttle position of sphere A 

- - - 

30 

ETAL(l) 

Dummy word 

- - - 

31 

ZETAL(l) 

Dummy word 

... 

2 to 56 

- - - 

Repeats words 7 to 31 for sphere B 

- - - 
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TABLE H.- OUTPUT FROM ASTAC MODEL 

Output during each engagement 

1. Initial position and velocity of each aircraft 

2. Position, velocity, thrust, drag, load factor, and attitude of each aircraft, and 

range and closure rate between aircraft, at preselected time intervals 

3. Time each new tactical maneuver is selected, and type of maneuver selected 

Output after each engagement (for each aircraft) 

1. Percent of time at - 

load factor: less than 2.5, 2.5 to 5.0, over 5.0 
Mach: less than 0.8, 0.8 to 1.2, over 1.2 

altitude: less than 10 000, 10 000 to 20 000, 20 000 to 30 000, over 30 000 feet 

2. Time in offensive position 

3. Fuel remaining 

4. Type and number of missiles carried and fired 

5. Cumulative probability of kill for gun, IR missiles, and radar missiles 

6. Time and distance of closest approach for each missile fired 

7. Time opponent was detected 

Summary for set of engagements (for each aircraft) 

1. Mean probability of kill for gun fire, IR missiles, and radar missiles 

2. Mean of altitude versus Mach versus load factor distribution, over all duels 

3. Mean percent of total gunfire time at 

load factor: less than 2, 2 to 4, 4 to 6, over 6 
range: 0 to 1000, 1000 to 2000, over 2000 feet 
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TABLE m.- QUESTIONS TO DETERMINE CELL OF SITUATION MATRIX IN AML 

(1) Is the opponent within a certain cone behind me? 

(2) Am I within a certain cone behind my opponent? 

(3) Is the opponent in an attitude and position so that he can fire at me ? 

(4) Am I in an attitude and position so that I can fire at my opponent ? 

(5) What is the rate of closure (closing fast, closing slowly, separating slowly, 
separating fast) ? 

(6) Is the line-of-sight angle (angle between the velocity vector and the line-of-sight 
vector) from the opponent to me less than 60° ? 

(7) Is the line-of-sight angle from me to the opponent greater than 90° ? 

(8) Is the rate of change of the line-of-sight angle from me to the opponent negative? 

(9) Is my specific energy rate P s greater than a given constant? 
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TABLE IV.- OUTPUT FROM AML PROGRAM 
Output during each engagement 

1. Position, velocity, lift, drag, thrust, load factor, angle of attack, specific excess 

power, throttle position, and attitude of each aircraft, and range and closure 
rate at selected time intervals 

2. Currently accumulated offensive, defensive, and gun envelope time 

3. Current AML cell state 

4. Trial maneuver data including number of maneuvers tested, load (g units) pulled 

during maneuvers, P s values for each maneuver, cell state for each maneuver 

Output after engagement 

1. Time on offense and defense for each aircraft 

2. Time in gun envelope for both aircraft 

3. Percent of run time on offense and in gun envelope 

4. Percent of offensive time in gun envelope 

5. Total run time 

6. Plots of aircraft trajectories projected on X-Y plane 

7. 3-D plots of each run (perspective and projected plots) 
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TABLE VI.- TIME LISTING OF SPECIFIC AIRCRAFT INFORMATION OUTPUT BY TLIST (TYPE II OUTPUT) 
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TABLE VII.- TIME LISTINGS OF ANGULAR RELATIONS BETWEEN TWO AIRCRAFT OUTPUT BY TLIST (TYPE m OUTPUT) 
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TABLE VIII.- SPECIAL EVENT MARKERS OUTPUT BY SUBROUTINE FLAGER 
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198.0 


TABLE IX.- HISTOGRAMS GENERATED BY TTHIS 


Range - Lambda A 

Lambda A (Lambda B) - Range 

Lambda A (+Range rate) - Range 

Lambda A (-Range rate) - Range 

Range - Altitude A 

Lambda A - Altitude A 

Lambda A - Altitude B 

Altitude A - Altitude B 

Velocity A - Altitude B 

Velocity A - Altitude A 

Normal Load A - Altitude A 

Normal Load Factor A (Mach A) - Altitude A 

Range - Lambda B 

Lambda B (Lambda A) - Range 

Lambda B (+Range rate) - Range 

Lambda B (-Range rate) - Range 

Range - Altitude B 

Lambda B - Altitude B 

Lambda B - Altitude A 

Altitude B - Altitude A 

Velocity B - Altitude A 

Velocity B - Altitude B 

Normal Load B - Altitude B 

Normal Load Factor B (Mach B) - Altitude B 



TABLE X.- THREE-WAY HISTOGRAM (X A , X B , R) OUTPUT BY SUBROUTINE TTHIS 
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TABLE Xffl.- FIRST PAGE OF S1RUN OUTPUT 
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TABLE XIV.- SECOND PAGE OF S1RUN OUTPUT 
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P PAL Lf.AD 0 TC 2.5 2.5 TO 5 5... IGRCS) 



TABLE XVI.- SI RUN GUNFIRE AND MISSILE LAUNCH OUTPUT 
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PLAN E IN SPHEKF A LAUNCHES MISSILE TYPI f IN ENVELOPE AI TIME 


TABLE XVH.- TYPICAL FIRST MOMENT TIME LISTING OUTPUT BY M123S 
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TABLE XXIV - TYPICAL PROGRAM MULTE OUTPUT GIVING MEANS OF IMPORTANT PARAMETERS 
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SECS IN MISSILE A ENVELOPE .5 .8 1.5 2.1 7.5 9.7 15.5 39.2 65. A 85.5 123.8 164 
SECS IN MISSILE B ENVELOPE .4 .* .8 1.2 . 4.7 6.1 13.8 33.9 58.6 77.9 113.0 150 
SECS IN MISSILE C ENVELOPE 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0.0 0 



TABLE XXV.- TYPICAL STANDARD DEVIATION OUTPUT BY PROGRAM MULTI2 
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TABLE XXVI.- TYPICAL t-TEST OUTPUT BY PROGRAM MULTI2 
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TABLE XXVIII.- TYPICAL THREE-WAY HISTOGRAM OF A*, X B , AND R OUTPUT BY PROGRAM MULTI2 
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TABLE XXX.- TYPICAL MULTIRUN OUTPUT FROM PROGRAM BZONE 
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Langley differential maneuvering simulator (DMS). 














Performance Data 



Figure 3.- DMS specific excess power programs. 
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Figure 6.- Typical output of program EXTRA. 
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Multirun Program 
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Figure 9.- Programs using DMS results. 
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Figure 10.- Main subroutines used in HASTE and types of output provided. 














Figure 12.- Time histories of range and off-boresight angle for both 
aircraft A and B provided by subroutine GRAFIC . 


RlTITuCSE Of SPHERE R 
UOOOS 0* FEET ) 




SPHERE fl: PILOT H -AIRCRAFT 6 / SPHERE B: PILOT 13 -AIRCRAFT 2 
DATA EVERY 10-0 SECONDS FROM 0-0 TO 233.0 SECONDS 

R OFFENSE. X DEFENSE. N NEUTRRL 





0-0 -2 »H -6 .8 1-0 1.2 1.4 1 . 

RUN 32 MRCH OF SPF€RE R 3-23-7? 

Figure 13.- Typical Mach-altitude plot output by program GRAFIC 











Figure 14.- Typical graphical representation of NADC missile program 
launch acceptance region (LAR) table. 
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Figure 15.- Typical Program 3-D output. 
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Figure 16.- Program BZONE zone definition. 
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